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In magnetized proto-neutron stars, neutrino cross sections depend asymmetrically on the neutrino 
momenta due to parity violation. However, these asymmetric opacities do not induce any asymmetric 
flux in the bulk interior of the star where neutrinos are nearly in thermal equilibrium. Consequently, 
parity violation in neutrino absorption and scattering can only give rise to asymmetric neutrino flux 
above the neutrino-matter decoupling layer. The kick velocity is substantially reduced from previous 
estimates, requiring a dipole field B ~ 10^^ G to get Wkick of order a few hundred km s~^. 
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Recent analyses of pulsar proper motion [Q indicate 
that neutron stars receive large kick velocities at birth 
{v]iick = 200 — 500 km on average, with possibly a 
significant population having velocity ^ 1000 km s~^). 
Evidence for high velocity neutron stars has also come 
from observations of pulsar bow shocks |^ and studies 
of pulsar-supernova remnant associations Support 
for natal kicks has come from the detection of geodetic 
precession m the binary pulsar PSR 1913-^16 ji), and 
' orbital precession in the PSR J0045-7319 binary and its 
fast orbital decay Q . In addition, evolutionary studies of 
the neutron star binary population imply the existence 
of pulsar kicks , and observations of nearby supernovae 
and supernova remnants also support the notion that su- 
pernova explosions are not spherically symmetric. 

Mechanisms for the pulsar kicks may be divided into 
two classes. The first class relies on hydrodynamical in- 
stabilities in the proto-neutron star |^ and nonspherical 
I perturbations in the precollapse core The asymme- 
tries in the density and temperature distributions nat- 
urally lead to asymmetric matter ejection and/or asym- 
metric neutrino emission. In this paper, we are concerned 
with the second class of models, where the kicks arise 
from asymmetric neutrino emission induced by strong 
magnetic fields. The fractional asymmetry a in the ra- 
diated neutrino energy required is a = iWwkickc/-E'tot 
{— 0.028 for Wkick = 1000 km s~^, neutron star mass 
M = 1.4 Mq and total neutrino energy radiated E'tot = 
3 X 10^3 erg) 

A number of authors have noted that parity viola- 
tion in weak interactions can lead to asymmetric neu- 
trino emission from proto-neutron stars [p|-pd|. It has 
recently been suggested [ p^ that the asymmetry in neu- 
trino emission may be enhanced due to multiple scatter- 
ings of neutrinos by nucleons which are slightly polarized 
by the magnetic field. Initial neutrino cooling calcula- 
tions in magnetic fields appeared to indicate that 
a dipole field of order lO^'' G is needed to produce kick 



velocity of 200 km s-\ These resuhs are wrong. 

Indeed, as we show in this paper, although the neutrino 
scattering cross-section is asymmetric, detailed balance 
requires that there be no cumulative effect associated 
with multiple scatterings in the bulk interior of the star 
where thermal equilibrium is maintained to a good ap- 
proximation. Moreover, we derive an explicit expression 
for the neutrino flux in a magnetic field. In addition to 
the usual diffusive flux, there is a drift flux (along the 
magnetic field) which is proportional to the deviation of 
neutrino distribution from equilibrium. Hence parity vi- 
olation can only induce an asymmetric neutrino flux near 
the surface of the star. 

Our starting point is the Boltzmann transport equa- 
tion for neutrinos (of a given species): 



at 



ft ■ V/,(k) 



ra/.(k)i 




[5/.(k)1 






+ 

SC 




abs 



(1) 



where k = fcJl is the neutrino momentum ($7 is a unit 
vector), /iy(k) is the neutrino distribution function (the 
position and time dependence are suppressed), and the 
scattering and absorption/emission collision terms are in- 
cluded on the right-hand side of the equation (we set 
Ti = c = ks — ^ throughout the paper). In this paper we 
concentrate on neutrino scattering by nucleons [v + N 
V -\- N) and electron neutrino absorption/emission from 
the processes n + ly^ ^ P + £ , P + i^e ^ n + . Other 
processes can be similarly considered, but are less impor- 
tant. 

First we examine the scattering rate: 
d^k' d^p d^p' 
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where W^^") = {2Tr)W^{P + K - P' - K) \M,,,{n,n')\^ 
is S-matrix squared divided by time, P and K {P' and 
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K') are the 4-niomenta of the initial (final) nucleon and 
neutrino, fl = /i/(k') is the final state neutrino distri- 
bution function and /jv and f'pj are the initial and final 
state nucleon distribution functions, respectively. Note 
that the energy of (nonrelativistic) nucleons includes the 
spin energy, —sfimB (where s — ±1 is the spin, fim is the 
magnetic moment) jTsd . In writing down eq. (^) we have 

implicitly used the relation W^/f^| = wj^^^j , reflecting the 
time reversal symmetry of the weak interaction. 

In the bulk interior of the proto-neutron star, neutri- 
nos are to a good approximation in thermal equilibrium, 
with fi, ~ Z^"-' (the Fermi-Dirac function). Using the 
explicit forms for the neutrino and nucleon distribution 
functions, together with energy conservation, we can ver- 
ify the equality 

(l-/("))(l-/Ar)/;./f ) =/(")/Ar(l-/;.)(l-/:(")), (3) 

representing detailed balance in thermal equilibrium. 
Therefore the only nonzero contribution to {df^/dt)sc 
must be proportional to the deviation from thermal equi- 
librium; there can be no drift flux along the magnetic field 
proportional to flf'^ . 

We now proceed to derive an expression for the neu- 
trino flux. Using the Fermi-Dirac distribution functions 
of nucleons, we can write the scattering rate (|^) in the 
form 
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where qq = k~k' is the energy transfer, dT/dk'dfl' is the 
differential cross-section (per unit volume) for scattering, 
as given by 
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and the "nucleon response function" is defined to be 
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where q= |k — k'|. Note that it is essential to retain the 
inelasticity in eq. (jj) in order to derive the neutrino drift 
flux. To lowest order in l/mjsf, the matrix element for 
iy~ N scattering jll, \Mss'{n, n')f is given by pi: 
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where we have used the interaction Hamiltonian as given 
in Ref. [nsl, G_f is the Fermi constant, ca and cy are 



the effective nucleon coupling constants, and A = ca/cv- 
Note that we have not summed over the initial or fi- 
nal state spins since both the nucleon distribution func- 
tions and the energy conservation delta function depend 
on spins. Time-reversal symmetry, |A/ss'(J^, J^')P ^ 
\Ms's{^', S^)P, can be explicitly verified for the matrix el- 
ement. The nucleon response function for zero magnetic 
field has been studied in Ref. |]l9| and can be directly gen- 
eralized to the B ^ case. Since the nucleon spin energy 
fj.mB is much smaller than other characteristic energies 
we may expand eq. (|^) to linear order in B. Combining 
the expressions for and Sss' with eq. (^, we find 
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where ji' = ft ■ O'. In (||), Aq corresponds to the B = 
result: 

Ao (k, k', fi') = A [(1 + 3A2) + (1 - A2) fj,'] 
I . / l + e-^° 
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and the corrections arising from nonzero B involve the 
coefficients 
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The reason for writing the cross section in the form of 
eq. (^) is that the angular dependence needed to find 
the moment equations is now manifest. This differential 
cross section shows the signs of parity violation through 
the • B and Vl' ■ B terms. Typically, SA± is smaller 
than ^0 by a factor of order fXmB/T. 

We now examine the macroscopic consequence of the 
asymmetric cross section in eq.(H). Let 
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where is the deviation of the spherically symmetric 
part of fn from the Fermi-Dirac value, while is the 
dipole dependent part which leads to the flux ||2l[. The 
first order moment equation is obtained by integrating 
eq. (0) against / dr2n/47r, with the result 
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The scattering term yields 



Here 
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leads to the usual B 
from the term 
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flux and the drift flux arises 



dk' j d^' [SA+{k,k',fi') 
+fi'SA_{k, k\ ^Ji')] [g,{k')C{k, k') + gAk)D{k, k')] , (17) 
where 
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Since gv(k) is negligible in the bulk interior of the star, 
5Vd (and hence the drift flux) is only important outside 
the neutrino decoupling sphere (outside of which there 
is no energy exchange between neutrons and matter). 
In this outer layer, the nucleons are nondegenerate and 
eqs. (^6|) and (^^ can be evaluated explicitly under the 
condition qq/T ^ 1. We find 
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where Kq'^'^^ is the zero-field scattering opacity (per unit 
volume), as given by 
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(n^ is nucleon number density), and the asymmetry pa- 
rameter is given by 
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Note that for nondegenerate nucleons qo/T ^ 
k/{mNTYl'^ <C 1. Naively one would expect from eq. (^) 
or eq. (|l7| ) that 5Yd would be suppressed by a factor of 
order q^/T However, inspection of eq. ( p^ ) shows 

that A± contains the term 2mNqvi/q^ which can be quite 
large and has the same parity (around k' = k) as qo/T. 
The result is that SY^ in eq. ( pi| ) is not suppressed by a 
factor of qo/T outside the neutrino decoupling sphere. 

The zeroth moment of eq. (^) describes the energy ex- 
change between matter and neutrino fields. The asym- 
metric scattering due to the magnetic field modifies this 



equation, but this only affects the rate at which the neu- 
trinos approach equilibrium. 

We now consider the effect of asymmetric absorp- 
tion/emission. For concreteness, we focus on the process 
Ve + n ^ p + . The absorption/emission rate can be 
written as 
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where s„, Sp are the neutron and proton spins, 
stands for summing over the electron phase space |2|] 
and l^(abs) ^ i^(abs) |2 j^E^-Ep- E^) is the tran- 
sition rate for absorption (with matrix element M'^'^^^^). 
Expanding /^(k) as in eq. (Oh, and using the equality 
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where the absorption opacity is given by 
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Calculation [ p^ , p7[ shows that can be expressed in 

terms of the B ~Q opacity ¥1^^"°^^ via 
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Taking the first moment of eq. (pj), we find 





[a/.(k)1 


/ 47r 





_ ^(abs) 



h.(fc) 



abs 



(29) 



This clearly indicates that in the bulk interior of the star, 
where g^ ~ 0, there is no drift flux associated with asym- 
metric absorption; the asymmetry from neutrino absorp- 
tion is exactly cancelled by that from emission p^ . 

Combining eqs. (|l|), (|l|), (|ol), (|ll), and (|§), we flnd 
the explicit expression for the neutrino flux in the outer 
layer of the star: 
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where k, 
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wiicic v*^; — ~^ '^\^^^\ and we have neglected 
dh.y/dt. Equation ( po| ) clearly shows that the neutrino 

drift flux (the terms with B) is proportional to 171,, which 
decreases rapidly with increasing optical depth. This in- 
dicates that asymmetric neutrino opacities due to parity 
violation do not affect neutrino transport in the bulk in- 
terior of the star. Asymmetry in neutrino flux can only 
arise from regions above the neutrino decoupling layer. 

Clearly, to determine the asymmetry in neutrino 
emission from a magnetized proto-neutron star requires 
knowledge of the temperature profile and the function 
gv{k) near the stellar surface. In principle, one can use 
the moment equations to calculate gu{k), but the temper- 
ature profile must also be determined self-consistently. A 
full neutrino transport calculation is beyond the scope 
of this paper. We now give a rough estimate of the 
neutrino asymmetry due to the i^e, Ve drift flux (the 
net drift flux associated with /i and r neutrinos is zero 
[p^ ). The asymmetry parameter for scattering is of or- 
der Csc O.OOBBis/T, where Bi^ is the field strength 
in units of lO'^^ G, and T the temperature in MeV. The 
asymmetry parameter for absorption, Cabs, is dominated 
for low energy neutrinos (^15 MeV) by the contribution 
from electrons in the ground Landau level |l0| , [l7|] 
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(where E,y is the neutrino energy in MeV, Q = 
1.29 MeV), and for high energy neutrinos by nucleon 
polarization (~ ^mB/T). The electron neutrinos de- 
couple from matter near the neutrinosphere, where typ- 
ical density and temperature are p ~ 10^^ g cm"'^, and 
T ~ 3 MeV. For a mean Ve energy of 10 MeV, eabs is 
greater than esc- The asymmetry in the Ve, Ve flux is 
approximately given by the ratio of the drift flux and 

the diffusive flux, of order eabs['*o'^^''''/'^o*°*'']- Averaging 
over all neutrino species, we find the total asymmetry in 
neutrino flux a ^ 0.2eabs- To generate a kick of a few 
hundreds per second would require a dipole field of order 
IQis G. 

After our paper was largely completed, a preprint by 
Kusenko et al. came to our attention, where the au- 
thors came to a similar conclusion. However, they only 
showed that there is no flux in equilibrium, and did not 
derive an expression for the flux. In addition, they con- 
tend that unitarity necessarily implies that the cross sec- 
tion cannot depend asymmetrically on the neutrino mo- 
menta, which is false (see eq. (||)). 
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